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Abstract 

We estimate the upper frequency cutoff of the galactic white dwarf binaries gravitational wave 
' background that will be observable by the LISA detector. This is done by including the modulation 

o ■ 

' of the gravitational wave signal due the motion of the detector around the Sun. We find this 

(N , 

. frequency cutoff to be equal to 10~^''^Hz, a factor of 2 smaller than the values previously derived. 
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pLn \ This implies an increase in the number of resolvable signals in the LISA band by a factor of about 

^' 4. 

^ ■ Our theoretical derivation is complemented by a numerical simulation, which shows that by using 

cn . 

the maximum likelihood estimation technique it is possible to accurately estimate the parameters 
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' of the resolvable signals and then remove them from the LISA data. 
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The gravitational waves from short-period binary systems containing white dwarfs and 
neutron stars are signals guaranteed to be observed by the Laser Interferometer Space An- 
tenna (LISA) mission. Recent surveys indicate that there exist about twenty of such systems 
that are emitting gravitational radiation of frequency falling into the LISA band. Popula- 
tion studies have also shown that the number of such sources will be so large to produce a 
stochastic background that will lie significantly above the LISA instrumental noise in the 
low-part of its frequency band. It has been shown in the literature (see for a recent study 
and for earlier investigations) that these sources will be dominated by detached white 
dwarf - white dwarf (wd-wd) binaries, with 1.1 x 10^ of such systems in our Galaxy. The 
detached wd-wd binaries evolve by gravitational-radiation reaction and the number of such 
sources rapidly decreases with increasing orbital frequency. It is therefore expected that 
above a certain limit frequency one will be able to resolve individual signals. 

In this communication we calculate this limit frequency by including the effect of the 
signal modulation induced by the motion of LISA around the Sun. Earlier derivations 
assumed that one would be able to resolve only one signal per frequency bin. However, the 
modulation of the signal depends on the position of the source in the sky and therefore, by 
using matched-filtering technique, it is possible to resolve signals with same frequency but 
incoming from different directions. 
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The number distribution, of the detached wd - wd binaries has been estimated in 
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where yr is expressed in seconds, and fo = 10"^'^ Hz. Thus, for an observation time of n 
years, the number of sources in one frequency bin,A^s, can be approximated by the following 
formula 
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If the signals were monochromatic one would expect to resolve roughly one source per 
frequency bin. Consequently, by equating A^^ to 1 and assuming 1 year of observation time, 
we would get an upper frequency cutoff equal to fo- However, by accounting in the data 
analysis for the modulation of the signal due to the motion of LISA around the Sun, it is 
possible to resolve signals of equal frequency but incoming from different sky positions. 
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To estimate this effect we represent the LISA response, to a monochromatic gravitational 
wave signal, in the following form 



h{t) = AoCos$(t) , (3) 

where the amplitude Ao is assumed to be constant, and the phase which includes the 
Doppler modulation, is given by 

^{t) =ujt+— cos(A) cos{nt - f3-(j)o) + $o. (4) 

Here uj is the angular frequency of the wave, Q = 27r/lyr, (A,/5) are ecliptic coordinates 
of the source, R = lAU, 0o is a known phase determining the position of the detector in 
the orbit around the Sun and $o is an unknown, constant phase of the signal. This is a 
reasonable approximation since the time scale over which the detector's response amplitude 
changes is significantly longer than that of the phase. 

It is convenient to introduce the following linear parameterization of the signal 

h{t) = Ao cos[$o + ujt + A cos{nt) + B sin(Ot)] , (5) 

where 

ujR 

A = cos A cos(/? — 0o), (6) 

iuR 

B = cos A sin(/5 — 0o). (7) 

In this parameterization the phase of the signal is a linear function of the signal's unknown 
parameters ($0; A, B). The reduced normalized correlation function of the signal is given 
by 3 

C{Aw, AA, AB) = ( f cos[Awx + AA cos(27rnx) + AB sin(27mx)] dxf + 

Jo 

( / sm[Awx + AA cos(27mx) + AB sin(27ma;)] dx)'^, (8) 
Jo 

where we have introduced a new dimensionless frequency parameter w = lyr n u, and the 
new dimensionless time variable x = t/(lyr n). The correlation function above is derived 
under the assumption that, over the bandwidth of the signal, the spectral density of the noise 
can be considered constant. As a consequence of the linearity of the phase of the signal in the 
parameters, the correlation function of two signals depends only on the difference between 



the values of the parameters and not on their absolute values. In the above parameterization 
the parameter space can be visualized as a truncated cone whose base and top are discs of 
radii and ^ respectively, with uo^ and uji being the angular frequencies of the upper 
and the lower edge of the bandwidth of the detector. For a given angular frequency uj the 
parameter space is a disc of radius Each point of the disc corresponds to two positions 
in the sky which differ by the sign in the ecliptic latitude angle A. 

In order to derive the background upper frequency cutoff we need to estimate the number 
of signals resolvable in a given frequency bin. This can be done by studying the correlation 
function between the signals in the A-B plane, and choosing a particular value for it. In this 
paper we have assumed two signals to be resolvable if their normalized correlation is less 
than 1/2. This choice of course is not optimal, and the final result of the upper frequency 
cutoff will depend on it. However, larger values of the correlation function imply a lower 
upper frequency cutoff, as it will be shown below. Here we do not identify the optimal value 
of the correlation function, which will be estimated in a forthcoming paper via Monte Carlo 
simulations. The reader should keep in mind this point when interpreting the conclusions 
of this paper. 

The correlation surface in the A-B plane can be approximated by an ellipse whose area, 
vaB) can be calculated using the covariance matrix (which is the inverse of the Fisher 
information matrix). The Fisher matrix, Fj^, for the signal (0) is given by 
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where z, j = ($o, w, A, B). Its expression has been derived under the assumption that, over 
the bandwidth of the signal, the spectral density of the noise can be considered constant. 
If we also assume the integration time to be a integer multiple of 1 year, then the Fisher 
matrix becomes 

/ 1 1/2 \ 

1/2 1/3 -l/(2n7r) 
1/2 

\ -l/(2n7r) 1/2 
In the above expressions for Fjj we have also normalized the signal-to-noise ratio to 1. 
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Since the area of the correlation elhpse is given by 

2 ^Tx'^n^ — 6 

where is the 2 by 2 sub- matrix of the covariance matrix for the parameters A and B, it 
follows that the number of resolved signals in a frequency bin centered on the frequency / 
is given by the ratio of the area of the disc of radius and the area of the correlation 

ellipse vab- In mathematical terms we have: 

Nr = No[j^ ' ^^^^ 

with 

■K{2-KfoR/cf , , 

iV„ = ^ — Jo I > _ (13) 

Vab 

If we equate the number of resolved signals, A^r, to the number Ng of expected signals in 
one frequency bin (Eq. Q), we get the following formula for the upper frequency cutoff, /,., 
of the background 

^ (Nouyn-'- ^^^^ 

Equation (fT^ implies a frequency cutoff fr = lO''^ '^ Hz when a year of observation time is 
assumed. 

By using the expression of the number of signals per frequency bin given in equation 
(^, we find that the number of detached white dwarf binaries above the frequency fr is 
around 6.8 x 10^. Although all these binaries can in principle be resolved, not all of them 
can be detected as some will have an amplitude smaller than the value of the detection 
threshold. To calculate the number of detectable binaries we first need to estimate the 
threshold corresponding to a satisfactory confidence of detection. We assume that we shall 
process the data by matched filtering. For the linear phase model of the signal given above 



we can adopt .nethods developed .1^. Ki.t we need to ea.enlate tKe n„.be. of eeUs, 

in the parameter space. These are the number of realizations of the optimal filter whose 
correlation is smaller than a pre-chosen value (1/2 in our case). This number is given by 
the ratio between the volume of the characteristic correlation hypersurface and the volume 
of the parameter space over which the search is performed. In our case the volume of the 
parameter space is given by 

TT /lAU^ ^ 



V = nYT-{ — ] U-uf), (15) 







(18) 



while the volume v of the correlation ellipsoid can be obtained by using the covariance 
matrix: 

71^2 det{C^) 

V = — . (16) 

Here C3 is the 3 by 3 sub-matrix of the covariance matrix for the parameters w, A, and B. 
In reference |4|, Eq.(77), it was shown that for Gaussian noise, the false alarm probability 
can be written as follows 

Pf = 1 - [1 - e-^«]^^ , (17) 
where To is the threshold on the following optimal statistics function P| 

I /.nyr 

n yr Sc 

In Eq. fITHjl Sc is the value of the two-sided power spectral density of the noise estimated in 
the center of the frequency bandwidth of the signal, and y(t) is the LISA data stream. 

Equation (fT7j) implies that, with a false alarm probability of 1 percent, the corresponding 
threshold signal-to-noise ratio is equal to 7.6. The corresponding number of resolvable white 
dwarf binaries with signal-to-noise above this threshold is equal to 3365, under the assump- 
tion of sources uniformly distributed in the galactic disk. As a comparison, the number of 
resolvable binaries calculated without including the effects of the phase modulation of the 
signals is equal to 919, a factor of 3.7 smaller than what we estimate. 

This threshold is of course rather conservative, since it reflects the assumption that we 
are searching for one or more signals in the data, and we should therefore regard it as an 
upper-bound. In reality these signals are present in the LISA data and, in order to calculate 
the lower bound for the threshold, we can assume each cell to contain a signal. Consequently, 
the false-alarm probability becomes 

Pf = e-^° , (19) 

since now = I in equation (fT7|) . The corresponding threshold signal-to-noise ratio goes 
down to 2.7, again for a false alarm probability of 1 percent. This implies that 47440 white- 
dwarf binaries of signal-to-noise ratio larger than this threshold will be resolvable. If we 
neglect the effects of the phase modulation of the signals, the number of resolvable binaries 
is equal to 13178 instead. 

As a demonstration of how well the parameters of the signals can be estimated by using 
the maximum likelihood method, we have performed a numerical simulation of our technique. 



We first estimate the parameters of the strongest signal, we remove it from the data, and 
then iterate these two steps until no signal crosses our predefined threshold. The numerical 
procedure for locating the maximum consists of two steps. The first is a coarse search 
of the maximum of the likelihood function on a predefined grid in the parameter space. 
This is then followed by a refinement around the region of the parameter space where the 
maximum identified by the coarse search is located. This second stage of the maximization 
is performed by using a numerical implementation of the Nelder-Mead algorithm, where the 
starting point of the maximization is determined by the values of the parameters identified 
by the coarse search. 

A graphical representation of the effectiveness of our method is presented in Figures 1 
and 2, which correspond to two difi'crcnt combinations of signals in the LISA data. Both 
figures show two sinusoidal signals of identical frequencies, and incoming from different two 
directions in the sky. The correlation between the two signals was chosen to be approximately 
equal to 1/2. In figure 1 we input signals with signal-to- noise ratios equal to 7 and 20, while 
in Figure 2 the signal-to-noise ratios have been increased to 20 and 60 respectively. The 
common frequency of the signals (/ = 1.5 x 10~^ Hz), and their incoming directions in the 
sky are unchanged in both figures. 

A visual inspection of the plots shows that the estimation of the parameters and the 
resolution of the signals (with the use of the maximum likelihood method) can be expected 
to be satisfactory for removing them from the data. This will make possible the search of 
other signals potentially present in the LISA data. 
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FIG. 1: Estimation and removal of one of two signals from the LISA data stream using a two-step 
maximum likelihood estimation. The two signals have signal-to-noise ratios equal to 7 and 20 
respectively, their frequencies are equal, but they originate from different points in the sky. The 
upper part of the figure shows the power spectra of the signals and the noise (dashed line), of the 
data after the stronger signal is removed (solid line), and of the noise only (dotted line). The lower 
part instead compares the shape of the signal injected (solid line) into the LISA data against that 
of the estimated one (dotted line). 
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FIG. 2: The signals are as in Figure la, but now the spectra shown include one in which the signals 
are both removed (solid line). 
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FIG. 3: As Figure la, but now the signal-to-noise ratios have been increased to 20 and 60 
tively 
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FIG. 4: The same as Figures lb, but the signals have signal-to-noise ratios equal to 20 and 60 
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